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The recently discovered FeAs-based superconductors show intriguing behavior and unusual dy- 
namics of electrons and holes which occupy the Fe d-orbitals and As 4s and 4p orbitals. Starting from 
the atomic limit, we carry out a strong coupling expansion to derive an effective hamiltonian that 
describes the electron and hole behavior. The hopping and the hybridization parameters between 
the Fe d and As s and p-orbitals are obtained by fitting the results of our density-functional-theory 
calculations to a tight-binding model with nearest-neighbor interactions and a minimal orbital ba- 
sis. We find that the effective hamiltonian, in the strong on-site Coulomb repulsion limit, operates 
on three distinct sub-spaces coupled through Hund's rule. The three sub-spaces describe different 
components (or subsystems): (a) one spanned by the d x 2_ y 2 Fe orbital; (b) one spanned by the 
degenerate atomic Fe orbitals d xz and d yz ; and (c) one spanned by the atomic Fe orbitals d xy and 
d z 2 . Each of these Hamiltonians is an extended t — t' — J — J' model and is characterized by different 
coupling constants and filling factors. For the case of the undoped material the second subspace 
alone prefers a ground state characterized by a spin-density-wave order similar to that observed 
in recent experimental studies, while the other two subspaces prefer an antiferromagnetic order. 
We argue that the observed spin-density-wave order minimizes the ground state energy of the total 
hamiltonian. 



PACS numbers: 74.70.-b,74.25.Ha,74.25.Jb,75.10.-b 

I. INTRODUCTION 

The recent observation of superconductivity in 
quaternary oxypnictidesi' 2 -' 3 -^'^'^, which are materials 
based on FeAs (and FeP) , has rekindled intense 
activity^^^^^^^^^^^^ 2 ^ 2 !^. to find a 
description of the strong electronic correlations present 
in these materials and in the cuprates which could be 
responsible for such phenomena. 

The structure of the new materials 2 (Fig.[T]) consists of 
FeAs layers sandwiched between LaO layers with rather 
weak interlayer coupling. The FeAs layers consist of 
a square lattice formed by the Fe atoms, while the As 
atoms, which sit just above and just below the plane (see 
Fig. [T|) , form FeAs4 octahedra squeezed along the c axis 
such that each Fe-As bond forms a ±30° angle with the 
Fe plane. 

As a function of temperature, the resistivity of the un- 
doped parent compound, which is not an insulator, shows 
a drop around 150 K 2 ^ before turning back up below 50 
K. In addition, the magnetic susceptibility also shows an 
anomaly at 150 K and it was argued 2 ^ that the parent 
material has a spin-density-wave (SDW) instability below 
150 K. Recent neutron diffraction studies 3 ^ demonstrate 
that the parent compound at 150 K undergoes a struc- 
tural distortion from tetragonal at high temperature to 
monoclinic at low temperature. Furthermore, these neu- 
tron diffraction studies show that below ~ 134 K, while in 
the monoclinic phase, it develops the SDW order shown 
in Fig. [T] Subsequent Mossbauer and muon spin rotation 
(fiSR) studies 3 -^ confirmed these findings: the structural 
transition is found to be around 156 K and the magnetic 



transition around 138 K. The structural distortion which 
brings about the monoclinic structure at low temperature 
is such that the rows of atoms which have their spins an- 
tiferromagnetically aligned are closer than the rows of 
atoms in the perpendicular direction. 

The electronic structure of LaO i- x F x FeAs has 
been studied by density functional theoryiii and 
by dynamical mean field theory^. There are argu- 
ments against phonon-mediated superconductivity^^ in 
LaOi-^F^FeAs. In addition, there are suggestions that 
a two-band mode l 16 i 17 ' 18 i 19 ' 22 may be the right effective 
hamiltonian to use in order to describe the low energy 
physics of these materials. 

Since the interlayer coupling is found to be weakii, 
in order to provide a simpler basis to understand the 
electronic structure of these materials we will focus on a 
single FeAs layer. As further justification of this choice, 
we have calculated the band structure of LaOFeAs and 
that of a single FeAs layer within density-functional the- 
ory and we find that the important features of the bands 
near the Fermi level obtained by the two calculations 
are essentially the same, including position relative to 
the Fermi level and overall dispersion features. Accord- 
ingly, our starting point is a hamiltonian which includes 
the five Fe <i-orbitals and four outer As orbitals (4s and 
Ap). In our interacting electron model we include direct 
Fe-Fe hopping, hybridization between the Fe d and the 
As 4s and Ap orbitals, the local Coulomb repulsion en- 
ergy for adding an electron on any of the Fe d states 
via Hubbard-type terms in the hamiltonian, and finally, 
coupling through Hund's rule. 

We determine the relative energy of the atomic orbitals 
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as well as the hopping and the hybridization matrix el- 
ements as follows: We carry out a set of first-principles 
calculations of the electronic structure of a single FeAs 
layer. Using as basis the above mentioned Fe and As 
states, we fit the results of the first-principles calculations 
using the tight-binding approximation to determine the 
values of the hopping matrix elements and the on-site en- 
ergy levels. We find that the energy difference between 
the atomic orbitals is less than 1 eV while the Coulomb 
repulsion to add two electrons on the same Fe d orbital is 
assumed to be significantly large r n i 12 ' 13 i 14 ' 15 i 28 . We per- 
form a strong coupling expansion in which the hopping 
and hybridization terms are used as perturbation and the 
unperturbed parts are the terms which correspond to the 
Coulomb repulsion energies for a pair of electrons placed 
on the same or different orbitals. Even in the case where 
these class of materials do not fulfil the requirement for 
a strong coupling expansion, the qualitative results or 
trends suggested by such a systematic analysis may still 
be useful. We systematically derive a low energy effective 
hamiltonian to describe this multi-band system. 

Through this analysis, we show that the relevant low 
energy degrees of freedom can be described by two 2- 
"flavor" subsystems, one in which the two flavors corre- 
spond to electrons in the d xy and d z i states of the Fe 
atoms and another in which the two flavors correspond 
to d xz and d yz Fe orbitals, and a third 1-flavor subsystem 
corresponding to the d x 2_ y 2 orbital. Each of those sub- 
systems is described by a t — t' — J — J 1 model where the 
spin degrees of freedom of all three subsystems couple 
through Hund's rule. Using this hamiltonian we are able 
to explain why the undoped material orders in the SDW 
pattern shown in Fig. [1] as reported by recent neutron 
scattering experiments^. 

The paper is organized as follows: Section II gives a 
detailed discussion of our density-functional-theory band 
structure calculations. Section III presents the strong- 
coupling limit analysis which leads to the effective spin- 
spin interaction hamiltonian derived in Section IV. Sec- 
tion V gives a discussion of the physics of the effective 
hamiltonian. 



II. ELECTRONIC STRUCTURE 
CALCULATIONS 

A. First- principles electronic structure 

Our first-principles calculations are performed within 
the framework of density functional theory (DFT) and 
the local density approximation (LDA) for exchange- 
correlation effects. We used the SIESTA code^ 2 - because 
it employs a localized basis of atomic-like orbitals for the 
expansion of the wavcfunctions which makes the interpre- 
tation of electronic wavcfunctions in the solid straight- 
forward and transparent, without the need for additional 
analysis such as projection to localized Wannier-type or- 
bitals. We use pseudopotentials of the Troullier-Martins 




FIG. 1: The structure of the FeAs layer in the FeAs based 
superconductors. The Fe atoms form a square lattice while 
the As atoms form two square sublattices one just above (or- 
ange) and the other just below (yellow) the Fe plane. The x 
and y axes used to characterize the orbitals are shown and 
the unit cell is indicated by the shaded square. The + and 
— signs on the Fe atoms denote the spin orientation observed 
in the neutron scattering experiment^. The labels are used 
in the text to explain the various types of antiferromagnetic 
exchange interactions. 



type^ to represent the interaction between valence elec- 
trons and ionic cores, and the Ceperley- Alder form for the 
exchange-correlation functional^. We have generated 
several different Fe pseudopotentials in order to check 
for any dependence of the results on this ingredient of 
the calculations. We find that the effect of the different 
pseudopotentials on the quantities reported in the fol- 
lowing, such as the electronic bands, is indiscernible. In 
addition to the SIESTA calculations, we have used the 
VASP code to ensure that there is reasonable agreement 
between the results of two very different computational 
schemes. The VASP code uses a plane-wave basis instead 
of localized orbitals^ and employs pseudopotentials of a 
different type^ than those in the SIESTA calculation. 

LaOFeAs, belonging to the tetragonal P4/nmm space 
group, has a layered structure 2 -. The FeAs layer serves 
as the carrier conduction channel and it has strong elec- 
tronic couplings within the layer. The unit cell for the 
simplified model system, a single FeAs layer, contains two 
Fe and two As atoms with a vacuum layer with thickness 
~ 19 A. The full system, bulk LaOFeAs has two atoms 
of each type (Fe, As, La and O) in the unit cell. 

In the SIESTA calculations, we choose an auxiliary real 
space grid equivalent to a plane-wave cutoff of 100 Ry, 
and use 8x8x1 Monkhorst-Pack k-point grid for the 
FeAs layer, 4 x 4 x 1 for the (2x2) FeAs supercell, and a 
8x8x4 grid for the LaOFeAs structure. For geometry 
optimization, a structure is considered fully relaxed when 
the magnitude of forces on all atoms is smaller than 0.04 
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FIG. 2: The band structure obtained from the LDA calculations, for the paramagnetic phase of: (a) the LaOFeAs compound; 
(b) the FeAs layer. In both cases the lines are from the SIESTA calculations and the points from the VASP calculation. The 
Fermi level is set at zero in each case, for the neutral material. 



eV/A. In the VASP calculation we use the same k-point 
grids and the default plane-wave cutoffs. 
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Extra Electrons per Fe 

FIG. 3: The total energy of the FeAs layer as a function of the 
added charge obtained with the VASP code. The reference 
energy is the energy of the neutral layer in the paramagnetic 
phase. 

The results obtained for bulk LaOFeAs and for the 
FeAs layer, using the SIESTA and VASP codes, for the 
paramagnetic phase are shown in Fig. O In both cases, 
the agreement between the two different computational 
schemes is remarkable. Moreover, there a very close sim- 
ilarity between the bands of the bulk LaOFeAs com- 
pound and the FeAs layer, especially in what concerns 
the features near the Fermi level. Notice that the neu- 
tral FeAs layer contains a different number of electrons 



per Fe atom than the bulk LaOFeAs compound. Ac- 
cordingly, we show the Fermi levels for the bulk and the 
layer slightly offset, to emphasize the similarity of the 
band-structure features. Moreover, by examining care- 
fully the wavefunction character of the bands near T and 
M, the points in the Brillouin Zone where electron and 
hole pockets appear in the bulk LaOFeAs compound, we 
establish that these states arise from orbitals associated 
with the Fe and As atoms. For these reasons, it is rea- 
sonable to concentrate on the band-structure of the FeAs 
layer alone, in order to build a comprehensive picture of 
the interacting electron system, presented in the follow- 
ing sections. 

We address next the issue of the spin configuration. 
We considered different spin configurations in the anti- 
ferromagnetic (AFM) phase. For these calculations we 
use a (2x2) FeAs supercell. The two spin configurations 
are: (i) AFM1, which is simply repetition of the spin 
configuration in the (lxl) unit cell; (ii) AFM2, which 
has the same spin alignment in one of the two diagonal 
directions of the Fe lattice and alternating spins along 
the other direction, as shown in Fig. [T] The total energy 
for the former spin configuration is lower by 0.95 eV per 
(2 x 2) cell, suggesting the AFM1 configuration is more 
stable than the AFM2 one for the neutral FeAs layer. 
We believe that this result is due to the fact that our 
calculation is for the charge neutral FeAs layer, where 
the Fermi level is lower than that of the LaOFeAs com- 
pound; namely, the FeAs layer as part of the LaOFeAs 
structure is negatively charged by an extra electron per 
Fe atom because the LaO layer is positively charged since 
the preferred oxidation state of La is La 3+ and that of 
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oxygen is O 2- . To show that this is the case, we carried 
out a calculation using the VASP code for a charged FeAs 
layer. In Fig.[3]the total energies of the AFM1 and AFM2 
phases are compared as a function of the added charge 
a per Fe atom. Notice that the AFM2 phase becomes 
energetically favorable for a ~ 1. The exact values of a 
can not be accurately determined in the context of these 
calculations because the energy differences are within the 
range of accuracy of DFT-LDA. The fact that the AFM2 
phase becomes the ground state for a = 1 will be estab- 
lished using the effective Hamiltonian derived in the next 
section, which captures the interacting-electron nature of 
the system in a more realistic manner, and as such, gives 
more reliable results for the magnetic phases. 



B. Tight-binding approximation model 
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FIG. 4: The band structure of FeAs layer in the paramagnetic 
phase obtained from: (a) the LDA calculations with the two 
different approaches (SIESTA - lines, and VASP - points); 
(b) the tight-binding approximation with nearest neighbor 
interactions and minimal orbital basis. 

Our goal next is to calculate a tight-binding hamilto- 
nian which approximately gives the same band structure 
as that of the first-principles results for the FeAs layer, 
using nearest-neighbor interactions only (hopping matrix 
elements) and a minimal orbital basis, consisting of the 

I 



We discuss the non-interacting parts, H a and T of the 
hamiltonian first. The operator dlfj^ creates an electron 



4s and Ap As orbitals and the 3d Fe orbitals. We expect 
that the hopping matrix elements, needed for the hamil- 
tonian upon which the strong coupling expansion will be 
based, are not significantly affected by the value of the 
filling factor (the position of the Fermi level). For the 
reasons discussed above, namely that a realistic picture 
of spin configurations can only arise from the interacting- 
electron treatment based on the effective hamiltonian, we 
will focus on reproducing with the tight-binding approx- 
imation the band structure of the paramagnetic phase, 
as obtained from the DFT-LDA calculations. 

Since introduction of electron doping is necessary in or- 
der to produce superconductivity in LaOFeAs-based ma- 
terials, we focus in reproducing as accurately as possible 
the features near and above the Fermi level. Note that 
hole-doping induced superconductivity in these materials 
has been also reported recently^. We also use informa- 
tion from the first-principles electronic wavefunctions to 
determine what is the optimal fit of the tight-binding 
approximation to the DFT-LDA results. As mentioned 
above, all the bands in the neighborhood of the Fermi 
level are associated with the As 4s and 4p and the Fe 
3d orbitals, and these features are well reproduced by 
the tight-binding approximation results. The best fit we 
could achieve is shown in Fig. [4j and compared to the 
first-principles results for an extended region near the 
Fermi level. The on-site and hopping matrix elements 
that produce this fit are presented in the next section, 
where these values are employed to construct the effec- 
tive hamiltonian. 



III. STRONG COUPLING LIMIT 

We consider the hamiltonian describing a single FeAs 
layer of Fe and As (or, in a more general formulation, P) 
atoms: 

H = H a +f + U, (1) 
with the three terms defined by: 

(2) 
(3) 

(4) 

I 

of spin a on the v th Fe d-orbital (y = 1,2,3,4,5 is the 



i.y.a l,a,<r 

t = - E + m - E E(cW + m> 

<ij>,cr,vu' i,a,i>,u 

u = E^E^kkuh E ^5E*kV) 

v i v.v'^v iaa' 

a i a.a'^a iaa' v,v' i 
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TABLE I: The on-site energies in eV for the Fe 3d orbitals as 
determined by approximating the results of the first-principles 
band-structure calculation, using the tight-binding approxi- 
mation discussed in Sec. [II] We also include the hopping ma- 
trix elements tw between two nearest neighbor Fe d-orbitals 
of the same type. The notation is explained in Sec. Mil 
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2 


3 


4 


5 


Fe 3d-orbital 


d x 2 — y 2 


dxz 


dyz 


dxy 


d z 2 




-4.6 


-4.5 


-4.5 


-4.5 


-4.2 


ti/V 


0.22 


0.5 


0.5 


0.43 


0.22 


TABLE II: Same as in Table Q] for the As 4s4p-orbitals. 
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2 


3 






As sp-orbital 


Px 


Py 


spt 








-4.7 


-4.7 


-5.2 





d T 2 — ? 



± \ - a\s) ±b\p z ). The LDA 
state together with the 



index that corresponds to the five d Fe orbitals 
d X z, d yz , d xy , and d z 2 respectively) at the site i which has 

an energy td{v)- The operator spf" creates an electron 
of spin a on the a th As which is one of three As orbitals. 
These As orbitals are formed as follows: first, because 
of the tetragonal symmetry, the two Ap x and Ap y orbitals 
remain unhybridized, while the 4s and Ap z As states form 
two linear combinations Ispf) 
calculation shows that the As sp 
Fe 4s form a bonding and an anti-bonding state, with the 
bonding state approximately 10 eV below the Fermi level 
and the anti-bonding state approximately 6 eV above the 
Fermi level. Therefore, these two states are not included 
in the tight-binding fit and the three As states included 
are the Ap x and Ap y orbitals and the sp+ hybrid. Thus, 
a = 1,2,3 corresponds to the cases of Ap x ,Ap y and sp+ 
respectively at the I th As site, with site energy e sp (a). 

V is the hybridization term between the Fe 3d orbitals 
and the As orbitals. The sum over means that it 
is over all four As sites I around the i th Fe site. The 
hybridization matrix element V^" is proportional to the 
wave function overlap of the v Fe d-orbital and the As 
a th sp-orbital. Some of these matrix elements are zero 
due to symmetry arguments and the most significant ones 
are of the order of, or less than, 1 eV, as obtained through 
our tight binding fit of the LDA results. 

In Tables IIIIII and in Table IIIII we give the non-zero 
matrix elements obtained by fitting the LDA results to 
the tight binding model (as explained in Sec. [II]) which 
includes the five Fe d-states and the three As 4s — Ap 
states for each of the two Fe and the two As atoms in the 
Fe2As2 unit cell, as well as the matrix elements t vv i and 
V va between these states. In addition to the above terms, 
the tight-binding approximation to the LDA results gives 
two hopping matrix elements t vv i\ the first for v — » d xz 
and v' — » d yz , which is t xzyx = 0.54 eV, and the second 
for v — * d xy and u 1 — > d z i, which is t xy z i — 0.20 eV. 
All other hopping matrix elements are either identically 
equal to zero due to symmetry or negligibly small. 



TABLE III: Same as in Table [I] for the hybridization matrix 
elements between Fe 3d- and As 4s4p-orbitals. The atom 
labels are those shown in Fig. [1] For the case of Fe atom 
labeled 2 in Fig. [T] the matrix elements are obtained from 
the same table by interchanging the labels 1 and 2 of the As 
orbitals and reversing the sign. 



yvc 


l(d x 2_ y 2) 


2(d xz ) 


3(d HZ ) 


A(d xy ) 


5(d z2 ) 


Px(l) 





0.1 





-0.2 





P. (2) 


0.4 


-1.45 








0.25 


P»(l) 


-0.4 





1.45 





0.25 


Pv(2) 








-0.1 


-0.2 





spf(l) 


-0.5 





0.7 





0.9 


spf(2) 


-0.5 


0.7 








-0.9 



We turn next to the interaction part, U, of the hamil- 



tonian. n v d (ia) 



^ier "2r 1S ^ ne number operator and 



Ul d) or Ui sp) give the Coulomb repulsion for a pair of 
electrons placed on the same d-orbital or the same s or 
p As-orbital. The llj^, (or U^) )is responsible for the 
Coulomb repulsion between different Fe-d (or As sp) or- 
bitals within the same atom. We will assume that the 
Coulomb-repulsion terms between the same or different 
Fe-d orbitals are significantly greater than their counter- 
parts for the As sp states, consistent with the general 
expectations for these values in the literature: is 
believed to be large of the order of 4-5 eV— iISiiiLHiSiSS 
while the parameter U^ sp ^ is expected to be much smaller 
than that. In addition, we will assume that the same site 
Fe-d Coulomb repulsion is larger than the inter-orbital 
Coulomb repulsion llj^j,. The term proportional to J H 
represents Hund's rule for the Fe d orbitals, with J H > 0, 
of order of less than 1 eV. We have neglected the Hund's 
rule coupling for As orbitals. 

Notice that the energy levels tdiy) and e sp (v) lie in 
the region —4.7 ± 0.5 eV, namely the energy difference 
between any pair of such states is less than 1 eV which 
is believed to be smaller than the characteristic Coulomb 
repulsion energy . In this paper we begin our analysis 
from the atomic or strong coupling limit, which implies 
that we have assumed that the energy scale U* defined 
below in Eq. [T7] is significantly larger than the hopping 
and hybridization matrix elements. If this condition is 
not fulfilled for this class of materials, it may still be 
instructive to discuss the qualitative features which a 
strong coupling expansion yields. In this limit the un- 
perturbed part of the hamiltonian Hq is 



H = H a + U, 



(5) 



and the hopping part T, which includes the hybridiza- 
tion, plays the role of perturbation. As discussed pre- 
viously, the FeAs layer in the undoped LaOFeAs parent 
compound has an additional electron relative to the neu- 
tral FeAs layer. This is due to the fact that LaO layer is 
expected to be in combined 1+ oxidation state. There- 
fore, the eight states considered above (five Fe d states 
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plus three As sp states) are occupied by 12 electrons. 
The atomic configuration is shown in Fig. [5] Since we 
have five more electrons than levels, five energy levels 
must be doubly occupied. Double occupancy of the Fe 
d orbitals is much more costly compared to the As sp 
orbitals; therefore, all three of the As orbitals must be 
doubly occupied and the only doubly occupied Fe d or- 
bital is d x 2_ y i, which has the lowest energy. The other 
four Fe d orbitals are singly occupied and the spin of these 
electrons are parallel because of Hund's rule, represented 
by the coupling Jjj . 
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FIG. 5: The occupation of the Fe and As atomic levels in the 
FeAs plane of the undoped parent compound. The ordering of 
the levels is shown schematically, as obtained from the tight- 
binding approximation parameters. 



Next we consider the effective hamiltonian which, in 
strong coupling perturbation theory, is given by 



H eff = EoPo+PoTfl, 

n = P + R(f &-&?&), 

6 

R 



Hex — En 



(6) 
(7) 

(8) 



where E$ is the ground state energy given by the energy 
of the state depicted in Fig. The H e ff operates in the 
subspace So formed by the degenerate ground states of 
Hq, that is, the subspace of states produced by the direct 
product of atomic states like those in Fig. [5] in which the 
spins of the four electrons occupying the four d-orbitals, 
from one Fe atom to the next, point either all up or all 
down. The operator Pq is a projection operator which 



projects into the subspace So and Q = 1 — P , that is, 
the operator which projects outside the subspace So- The 
above equation can be formally solved iteratively to yield 
the Raylcigh-Schrodinger expansion as a power series in 
T . The leading term is T which, when restricted in this 
subspace So, becomes just the direct Fe-Fe hopping terms 



PofPo = - 



(t >d {v '^d {v) 



h.c), 



(9) 



In the right-hand-side of the above equation we have 
omitted the projection operators by assuming that we 
will restrict ourselves to the subspace So- 
la the following, in order to simplify the calculation, we 
will take lli — U to be independent of v and we will as- 
sume that U is much larger than the atomic energy level 
difference (which was found to be less than 1 eV within 
our TB approximation) and significantly larger than the 
hopping and hybridization parameters. In addition, we 
will take Ui% = U, U { v sp) = U sp and U { a % = U sp , i.e., 
to be independent of v, v' (or a, a'). 



IV. EFFECTIVE SPIN-SPIN INTERACTION 
HAMILTONIAN 

A. Interaction between same-type Fe orbitals 

First, there are the familiar second order processes aris- 
ing from the direct Fe-Fe hopping through the matrix 
elements t uv which give rise to an antiferromagnetic ex- 
change interaction of the form 



rtv v 



7(2) _ 

<J i, u 



7(2) \^ qu . 3» 



At 2 

~1T' 



(10) 



(ii) 



In addition, we have antiferromagnetic coupling of the 
spins of two different type Fe <i-orbitals due to the hop- 
ping terms t XZ}VZ and t xy z 2 which give rise to 



fj jr(2) \ qU qv' 



T (2) _ ^ly 



(12) 
(13) 



where v, v' can be either xz, yz or xy, z 2 . These processes 
take place only between nearest neighbors such as the Fe 
atoms 1 and 2 in Fig. [TJ These second order contribu- 
tions are obtained from the square of the matrix elements 
listed in Table HVl by multiplying them with A/U. There 
are no second-order nnn contributions to the spin-spin 
interaction. 

The next terms to leading order, beyond the first and 
second order terms discussed above, are fourth order pro- 
cesses involving an Fe <i-orbital and the As sp-orbitals 
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TABLE IV: The second-order nn contribution to the spin-spin 
couplings is obtained from the matrix elements listed below 
(in units of eV 2 ) by multiplying them with A/U . 



t 2 , 

b vv' 


dx 2 — y 2 


dxz 




dxy 


d z 2 


d x 2 —y 2 


0.048 














dxz 





0.25 


0.29 








dyz 





0.29 


0.25 








dxy 











0.185 


0.04 


d z 2 











0.04 


0.048 



with which it hybridizes. The contribution of all sp or- 
bitals of the two possible intervening As atoms to the 
exchange interaction between the same d orbital of two 
nn Fe atoms is given by 3 -^ 



7(4) 



ra 



\ J" 

lu* 



2b 2 



U* 



b = 



U* + e 



where 



U* = U + W- (iU sp + U sp ), 



(14) 

(15) 
(16) 

(17) 



and U* is assumed significantly larger that all other en- 
ergy scales in the problem. Here 



,(a) - e d (v), 



(18) 



is the energy difference between the As sp state and that 
of the Fe d orbital. The subscripts 1 and 2 in the matrix 
elements refer to the fact that the two nn atoms 1 and 2 
are at 90° angle relative to the position of the intervening 
As atom and, therefore, V\ is from Tabic IIII1 while Vi 
from Table IIIII by interchanging the index 1 and 2 and 
the sign of the matrix elements. 

Since e < 1 eV for any combination of Fe-d and the 
three As sp orbitals and U* ~ 5 eV, if e is neglected in 
the above expression the error in our estimate for the ex- 
change couplings J will be rather small. We are going to 
use the full expression given by Eq . [Tol when we compute 
the coupling constants J between the same-type Fe or- 
bitals. We will also use this approximation of neglecting e 
because it simplifies the results and this gives additional 
insight. With this approximation, we find that 



7(4) 

'J 1/ n 



4^U 

u* 3 ' 

6 



(19) 
(20) 



a=l 



Since the value of U* is not known, we can use the ex- 
pression given by Eq. [201 as a measure of the relative 



spin-exchange interaction coupling. Namely, to obtain 
the actual values of Jut , and the J^J, with v ^ v' to 
be discussed next, we just need to multiply the values 
given in Table IVl by 4/(C/*) 3 . Therefore, for two nn Fc 
atoms, such as Fe atoms 1 and 2 in Fig. [T] the spin-spin 
interaction coupling constant J vv is obtained by adding 
the second and fourth order contributions jjh) and 
respectively. 

On the other hand, for two next nn Fe atoms, such as 1 
and 3 in Fig. [TJ there are no second order processes since 
there are no direct hopping between such atoms. The 
fourth order superexchange contributions is obtained as 
follows: 



J, 



la 
vv 



j! o(") ot") 

<<ij>> 

1 6 

- V j' a 

a=i 

W* U* + e 



(21) 
(22) 

b' = 2 ;7 ' . (23) 
U* + e v ' 



Notice, that in this case the same matrix elements V\ 
are involved because both atoms participating in the su- 
perexchange form the same angle with the intervening 
As atoms. In addition, there is a factor of 2 difference 
between the above expression and Eq. [16] because there 
is only one possible intervening As atom for two fixed Fe 
atoms. The simplified expression, when e is neglected is 
the following: 



7'(4) 

A 1 



U 



,3 ' 



^E(^r) 4 - 

a = l 



(24) 
(25) 



The values of the constant A' vv are given as the diagonal 
matrix elements of Table IVD 



B. Interaction between different-type Fe orbitals 

There is an effective spin-spin interaction between cer- 
tain Fe orbitals of different type, J vv > with v ^ v' . This 
type of nn and nnn spin-spin interaction for most of the 
orbitals is significantly smaller than the J vv coupling con- 
stants between the same orbitals. To obtain an estimate 
of these we use the simplified expressions where e is ne- 
glected. Namely, 



j(4) 



r'(4) 



4A„ 



U* 

6 

E 

a=l 
AA' 



(vrvr ) 2 , 



u 



* 3 



(26) 
(27) 
(28) 
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TABLE V: Fourth-order sum of matrix elements contributing 
to the nn spin-spin couplings. 





d x 2 — y 2 


d xz 


dyz 


dxy 


d z 2 


dx 2 — y 2 


0.13 


0.12 


0.12 


0.01 


0.40 


dxz 


0.12 


0.04 


0.24 


0.08 


0.40 


dyz 


0.12 


0.24 


0.04 


0.08 


0.40 


dxy 


0.01 


0.08 


0.08 





0.01 


d z 2 


0.40 


0.40 


0.40 


0.01 


1.31 



TABLE VII: Second-order terms contributing to the effective 
nn hopping. 



B vv i 


d x 2_ y 2 


dxz 


dyz 


dxy 


d z 2 


d X 2—y2 


-0.50 


0.31 


0.31 








dxz 


0.31 


0.29 


-0.49 


-0.29 


-0.65 


dyz 


0.31 


-0.49 


0.29 


0.29 


0.65 


d x y 





-0.29 


0.29 





0.1 


d z 2 





-0.65 


0.65 


0.01 


1.62 



TABLE VI: Fourth-order sum of matrix elements contributing 
to the nnn spin-spin couplings. 



A' 


dx 2 —y 2 


dxz 


dyz 


dxy 


d z 2 


d x 2 —y 2 


0.18 


0.46 


0.46 





0.42 


dxz 


0.46 


4.66 








0.53 


dyz 


0.46 





4.66 





0.53 


dxy 

















d z 2 


0.42 


0.53 


0.53 





1.32 



TABLE VIII: Second-order terms contributing to the effective 
nnn hopping. 





d x 2_ y 2 


dxz 


dyz 


dxy 


d z 2 


d x 2 — y 2 


0.82 


-0.93 


-0.93 








dxz 


-0.93 


2.60 





-0.02 


-0.99 


dyz 


-0.93 





2.60 


0.02 


0.99 


dxy 





-0.02 


0.02 


0.08 





d z 2 





-0.99 


0.99 





1.74 



1 6 

4-' = -E^r^r; 



(29) 



a=l 



The values of the constants A vv i and A' vv , are given as 



the off-diagonal matrix elements of Tables [V] and [VT 



C. Hopping between same-type Fe orbitals 

There are contributions to the effective hopping matrix 
elements due to second order processes. Namely, pro- 
cesses in which an electron from a doubly occupied As 
orbital momentarily hops to the nn Fe d orbital and then 
an electron from the doubly occupied doped d-orbital 
hops to the singly occupied As orbital left behind. These 
processes give rise to the following expression 



st v , v = 2J 



yavyav 

U* + e 



(30) 



where e is given by Eq. [THJ Again, since the value of U* 
is not known, for large enough values of U* (e < leV) 
we can neglect e in the above expression to obtain the 
following expression 



0Zv < v ~ jj* ' 



(31) 
(32) 



and the values of B ul/ are given as the diagonal elements 
in in Table IVIII The actual estimates for St uv can be 
obtained by dividing the values in the table by U* . The 
total effective nn hopping is given as 



(33) 



In the case of next nn such as the Fe atoms 1 and 3 in 
Fig Q] we obtain 



p D vv 

VU jj* ■ 



(34) 
(35) 



The values of B' vv are given as the diagonal elements of 
TableEEEH 



D. Hopping between different-type Fe orbitals 

There is a second order process by means of which a 
doubly occupied site can effectively hop to a nn Fe d 
orbital of different type by involving an intervening As 
sp orbital. These contributions are smaller than those 
connecting two of the same-type Fe d orbitals and they 
can be approximated by 



St. 



u* 

6 



(36) 
(37) 



The values of B„„/ are given by the off-diagonal elements 
of Table ECU 

There is also a second order process which gives rise to 
hopping between nnn Fe d orbitals of different type. For 
this case we obtain 



St' 



B 'y V i 

u* 

6 



(38) 
(39) 
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TABLE IX: The estimated matrix elements for the J vv i and 
J' vv , for U = 5 eV. 



r 

J 


d x 2_y2 


A 

&XZ 


A 

yz 


a 

u>xy 


rl i 

a z 2 


P 
j 












d x 2_ y 2 


n n4 


n n 

u.u 


n n 

u.u 


n n 

u.u 


U.Ul 




n 

u 


n 

u 


n ni 

U.Ul 


n m 

U.Ul 


U.Ul 


A 

u>xz 


n n 


u.zu 


n 94 


n n 

u.u 


n ni 

U.Ul 




0.01 


0.07 


o 


0.0 


0.01 


dyz 


0.0 


0.24 


0.20 


0.0 


0.01 




0.01 





0.07 


0.0 


0.01 


dxy 


0.0 


0.0 


0.0 


0.15 


0.03 



















d z 2 


0.01 


0.01 


0.01 


0.03 


0.08 




0.01 


0.01 


0.01 


0.0 


0.02 



The values of B' vv , are given by the off-diagonal elements 

of Table Em! 

The total effective nn and nnn hopping for v ^ v' are 
given as 

iuu> = W +St vv >, (40) 
~tv = Kw (41) 

E. Effective Hamiltonian 

In summary the effective hamiltonian is given by 

where each of the "Hyy terms above may be written as 

7~lv,v> = ~~ ^ ] tv,v'Cj u i a Civcr ~ y ] ^u,v' C ju' a C ivo 
<ij>,<J «ij»,<? 

+ "£j„, u ,S»-S? + J2 J ly^-S^. (43) 

<ij> «ij>> 

Next, we will provide estimates of the coupling con- 
stants t' vv i, J v ^i and J' v v , involved in the above 
model based on the values of the parameters obtained 
from fitting the LDA results to the tight binding model. 
The matrix elements for nn hopping t vu i and spin-spin 
interaction J vv i as well as their counterparts for next nn 
interactions, that is, between sites diagonally across in 
the square lattice formed by the Fe atoms, t' vv , and J' vv , , 
are given in Tables HXIXi for U* = U = 5 eV and in 
Tables IXIIXIII using U* = U = 3 eV. A more simpli- 
fied model than the one given above is discussed in the 
following section. 

V. DISCUSSION 

First, by examining the Tables HXIXIXIIXlTl we notice 
that to a reasonable degree of approximation the follow- 
ing three subspaces couple with each other rather weakly: 



TABLE X: The estimated matrix elements for the t vv i and 
t' vv , for U = 5 eV. 















L 


dx^ — y^ 


A 

UXZ 


yz 


A 

u>xy 


rl n 

a z 2 


z 












d x 2_ y 2 


n 1 9 

U. 1^ 


n nfi 

U.UU 


n nfi 

U.UU 


n 

u 


n 

u 




n ns 


n no 

-u.uy 


n no 

-u.uy 


n 

u 


n 

u 


rl 


n nfi 

U.UU 


n ^p, 


n 44 


n nfi 

-U.UU 


ni'i 
-u. io 




-0.09 


0.26 


o 


o 


-0.1 


dyz 


0.06 


0.44 


0.56 


0.06 


0.13 




-0.09 





0.26 





0.1 


dxy 





-0.06 


0.06 


0.43 


0.22 













0.01 





d z 2 





-0.13 


0.13 


0.22 


0.54 







-0.1 


-0.1 





0.17 



TABLE XI: The estimated matrix elements for the J„„/ and 
J' uv , for U = 3 eV. 



J 


d x 2 —y 2 


dxz 


dyz 


dxy 


d z 2 


J' 












d x 2 — y 2 


0.08 


0.02 


0.02 





0.06 




0.01 


0.03 


0.03 





0.03 


dxz 


0.02 


0.34 


0.42 


0.01 


0.06 




0.03 


0.35 








0.04 


dyz 


0.02 


0.42 


0.34 


0.01 


0.06 




0.03 





0.35 





0.04 


dxy 





0.01 


0.01 


0.25 


0.05 



















d z 2 


0.06 


0.06 


0.06 


0.05 


0.26 




0.03 


0.04 


0.04 





0.10 



(a) one spanned by the d x 2_ y 2 Fe orbital; (b) one spanned 
by the degenerate atomic Fe orbitals d xz and d yz ; (c) one 
spanned by the atomic Fe orbitals d xy and d z 2 . Notice 
that the most significant off-diagonal matrix elements arc 
those which couple the d xz and the d yz orbital and those 
which couple the d xy to the d z 2 orbital. There are other 
smaller off-diagonal matrix elements which couple these 
subspaces weakly. These three subspaces, however, are 
much more strongly coupled through J H . 

We consider the undoped (LaOFeAs) case first. Be- 
cause six electrons should occupy the five Fe d orbitals 
(see Fig. [5]), the lowest energy subspace spanned by 
d x 2_ y 2 is occupied by two electrons, and also each of the 
other two subspaces is also occupied by two electrons. 
As mentioned above, these three subspaces are coupled 
mainly because J H ^ 0. Furthermore, the bands formed 
in any given subspace are intersected by the bands formed 
in the other two subspaces because their atomic energy 
difference is small compared to their bandwidth. 

There is a great degree of magnetic frustration, as 
noted in Refs. 20,21,25[, especially in the subspace 
spanned by d xz /d yz . In this subspace the next nn 
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TABLE XII: The estimated matrix elements for the t vv i and 
t' vv , for U = 3 eV. 



i 

p 

i 


d X 2—y2 




A 


j 

Uxy 




d x 2_ y 2 


n (Ti 


n in 


n in 


n 


n 




n 1 a 


-u.iu 


n i fi 


n 

u 


n 


A 

u>xz 


n 1 n 


n fin 


n 


n i n 

-u. 1U 


n 99 




-0.16 


0.43 


o 


o 


-0.17 




0.10 


0.38 


0.60 


0.10 


0.22 




-0.16 





0.43 





0.17 


dxy 





-0.10 


0.10 


0.43 


0.23 













0.01 





d z 2 





-0.22 


0.22 


0.23 


0.76 







-0.17 


0.17 





0.29 



Nearest neighbor coupling/2 
Next nearest neighbor coupling 




% 0.3 



U (eV) 

FIG. 6: The calculated nn J X z,xz/2 is compared to the nnn 
J'xz,xz as a function of U* (using U* = U). The condition 
for the instability of the (it, tt) order to the columnar order is 
J'xz,xg > Jxz,xz/2 which occurs for values of U* < 4.3 eV. 



coupling J' becomes greater than J/2 for any value 
of U* < A.ZeV (using U* = U)(See Fig. |BJ). When 
J < 2J', the observed columnar antifcrromagnetic or- 
dering is favorabl o 21 ' 25 relative to the familiar (tt, 7r) an- 
tiferromagnetic order. The subspace spanned by d x 2_ y 2 
has net spin zero and, therefore, is not expected to con- 
tribute significantly to the choice of magnetic order. On 
the other hand, there seems to be less frustration in the 
subspace spanned by d xy /d z 2, which is half-filled; we 
therefore expect that a long-range antifcrromagnetic or- 
der should characterize the ground state of this subspace 
if it were uncoupled from the d xz /d yz subspace. However, 
due to Hund's rule coupling J H , the spin orientation of 
all subspaces should be common. The conflicting prefer- 



ences of these two subspaces, which are forced to make a 
common choice, introduces further frustration of relative 
spin orientation. 

We expect that the subspace spanned by d xz /d yz drives 
the system to a global columnar orde r 20 ' 21 ! 25 because it 
is characterized by the larger couplings. The presence of 
a large J' in the subspace spanned by d xz /d yz might im- 
pose the observed columnar order through the relatively 
large Hund's rule coupling Jh ~ 0.5eV '. The fact that 
the subspace d xz /d yz prefers the columnar order and the 
subspace d z 2 jd xy prefers the (tt, tt) order creates frustra- 
tion which may also explain the fact that the observed 
moment per Fe atom is small. 

We emphasize that unlike the case of undoped 
cuprates, the undoped parent compound in the case of 
the oxypnictides is not an insulator. As can be inferred 
from Fig. [T] the motion along the ferromagnetic direc- 
tion is not hindered and, therefore, the undoped material 
is expected to demonstrate anisotropic transport in the 
SDW phase. 

We would like to discuss the case of the neutral FeAs 
layer which was considered in Sec. [Til where it was found 
that the ground state is characterized by (tt, tt) order. 
The case of the neutral layer has five electrons per Fe 
atom and this implies that all Fe d orbitals should be 
singly occupied. Therefore, the subspace spanned by 
d x 2_ y 2 is no longer characterized by spin zero. This 
means that the subspace d xz /d yz in order to drive the 
columnar order has to compete against not just one but 
two subspaces which prefer the (tt, tt) order. 

The phenomenological hamiltonian considered in 
0j|2ll| and in Ref. [25[ to introduce frustration, is 
different from the one we derived based on a more rig- 
orous approach, which is more complex. The next step 
would be to study the hamiltonian given in Eqs. P2"]) . p5]) 
by various analytical and numerical techniques, which is 
beyond the scope of the present work. 

While the estimated nn antifcrromagnetic coupling 
constants are of similar magnitude to the one in the 
cuprous oxides^, in the oxypnictide materials there is 
magnetic frustration mainly due to the fact that the nnn 
antiferromagnetic coupling for the d xz /d yz subspace is 
large. Therefore, assuming that the pairing interaction 
between electrons is of magnetic origin, it is not clear if 
the pairing energy scale is larger or smaller compared to 
that in the cuprate superconducting materials. The pair- 
ing energy scale in the present model may be enhanced 
by the "flavor" factor, that is, the number of states span- 
ning the subspace where the added electrons go in the 
case of electron doping, and by the fact that the hopping 
and spin-exchange matrix elements are estimated to be 
somewhat larger compared to those in the case of the 
cuprates (see Tables IDCIXIXIIXlTf . Therefore, it is con- 
ceivable that this new class of superconductors could lead 
to higher critical temperatures upon future optimization 
of the doping agents and other factors. 

A very important difference between the oxypnictides 
and the cuprates is that the five-fold sector can be 



_ Ref. 2 



11 



thought of as formed by three subsectors, two 2-flavor 
sectors and a third 1-flavor sector. The d xz /d yz sec- 
tor prefers the SDW order depicted in Fig. [TJ the sec- 
tor spanned by the d x 2_ y i orbital has spin zero and the 
other sector spanned by d xy /d z 2 prefers antiferromag- 
netic long-range order. These subspaces are coupled by 
Hund's rule which, we believe, leads to the SDW order 
with ferromagnetic order along one direction and anti- 
ferromagnetic ordering between such chains. As in the 
case of cuprates, superconductivity in the oxypnictide 
materials might coexist with SDW order— but these are 



expected to be to some extent competing orders as found 
in neutron^ and fiSR studies^ done on the supercon- 
ducting doped materials. 
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